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Abstract 

A criterion of thermal stability is derived for electrically charged quantum black holes having large horizon 
area (compared to Planck area), as an inequality between the mass of the black hole and its microcanon- 
ical entropy. The derivation is based on key results of Loop Quantum Gravity and equilibrium statistical 
mechanics of a grand canonical ensemble, with Gaussian fluctuations around an equilibrium thermal con- 
figuration assumed here to be a quantum isolated horizon. No aspect of classical black hole geometry is 
used to deduce the stability criterion. Since, no particular form of the mass function is used a priori, our 
stability criterion provides a platform to test the thermal stability of a black hole with a given mass func- 
tion. The mass functions of the two most familiar charged black hole solutions are tested as a fiducial check. 
We also discuss the validity of the saddle point approximation used to incorporate thermal fluctuations. 
, Moreover, the equilibrium Hawking temperature is shown to have an additional quantum correction over 

O" 1 the semiclassical value. 

1 Introduction 

The idea of a quantum black hole has been articulated longer than a decade ago within Loop Quantum Gravity 
[1], [2], [3], [4]. This is an effective description wherein a three dimensional SU(2) Chern Simons Theory [1], [5], 
[6] coupled to point-like sources on punctures made by edges of spin networks describing bulk quantum geometry, 
is shown to describe a quantum isolated horizon [7], [9] or, in other words, a null trapping horizon [11]. The 
dimensionality of the Hilbert space of the Chern Simons theory leads, for large horizon areas (in comparison to 
the Planck area) to the Bckenstein-Hawking area law together with universal quantum corrections [3]. These 
results have been rigorously reviewed and redcrived recently by a number of authors [12], [13], [14]. 

Quantum black holes not isolated from an ambient thermal reservoir have also been considered in the 
past [15], [16-18], [19]. In this approach one uses certain key results of Loop Quantum Gravity like the discrete 
spectrum of the area operator [27, 28] and the central assumption that the thermal equilibrium configuration is 
indeed an isolated horizon whose microcanonical entropy, including quantum spacetime fluctuations have already 
been computed via Loop Quantum Gravity. The idea here has been the study of the interplay of thermal and 
quantum fluctuations, and a criterion for thermal stability of such horizons has been obtained [18], [19,20], 
using a 'thermal holographic' description involving a canonical ensemble and incorporating Gaussian thermal 
fluctuations. The generalization to horizons carrying charge has also been attempted, using a grand canonical 
ensemble, even though a somewhat ad hoc mass spectrum has been assumed [17]. 

In this paper, we attempt to redcrive a thermal stability criterion for charged quantum horizons, without 
any ad hoc assumptions on the mass spectrum. With the benefit of hindsight, arguments which place the earlier 
formulation on a more solid footing are presented, together with novel aspects which enable us to sidestep earlier 
restrictions. A comparison with semiclassical thermal stability analyses of black holes [22] is made wherever 
possible. The range of validity of the saddle point approximation around the equilibrium configuration is 
examined to ensure the self-consistency of the Gaussian approximation. 

The paper is organized as follows: In Subsection(2.1) there is a short account on the mass associated with 
the horizon of a black hole. In the rest of Scction(2) we review the idea of Thermal Holography as laid out 
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in [20]; the primacy of the boundary partition function of a Grand Canonical ensemble, in situations where 
the bulk Hamiltonian is a constraint, is established. The boundary partition function is then evaluated in 
Section(3) within the saddle point approximation, choosing an isolated horizon as the equilibrium configuration. 
Valid existence of the saddle point is shown to lead to the thermal stability criteria in terms of a second 
order partial differential inequalities involving the equilibrium mass and the microcanonical entropy. Solving 
one of those inequalities yields a more accessible form of the stability criterion in terms of a competition 
between 'energy-driven' and 'entropy-driven' aspects of the system. In Section(4) the range of validity of the 
saddle point approximation for the different cases is discussed, thermal stability criterion with semiclassical 
thermal stability properties of Reissner- Nordstrom and Anti-de Sitter Reissner- Nordstrom black holes is given. 
Scction(5) contains our concluding remarks. 

2 Thermal Holography 

2.1 Horizon Energy 

For a consistent Hamiltonian evolution for spacetimcs admitting internal boundaries (isolated horizons, repre- 
senting black holes at equilibrium) there must be a first law associated with each internal boundary (6), assumed 
to be a null hypersurfacc with the properties of a 'one-way membrane' [7, 9] given by 

8E\ = tr5A b + &6Q b 

where E\ is the classical energy function associated with the horizon, k* and <!>' are the surface gravity and 
the electric potential respectively of the horizon, Q b is the horizon charge. All the quantities are defined for a 
particular choice of time evolution vector field t M . The family of time evolution vector fields [i^] satisfying such 
first laws on the horizon are the permissible time evolution vector fields. These evolution vector fields also need 
to satisfy other boundary conditions. Each of these time evolution vector fields associates a classical energy 
function with the horizon which is a function of area and charge for Einstein-Maxwell theory. In arbitrary 
non-stationary cases (radiation may be present arbitrarily close to the horizon) for a particular time evolution 
vector field (t), the Hamiltonian formulation yields 

H = E ADM - M b 

where H l = Hamiltonian associated with the spacetime region in between the black hole boundary(6) and 
the boundaries at infinity, Ml = the mass associated with the horizon (b), E ADM = the usual ADM energy 
associated with the spatial boundary at infinity for a permissible vector field t M . H* is the Hamiltonian of 
the covariant phase space, which is the space of various class of solutions of the Einstein equations admitting 
internal boundaries. For stationary spacetimes the global timelike Killing field (£ M ) is the time evolution vector 
field. On physical ground one can say that there is nothing between the internal boundary and the boundary at 
infinity for stationary spacetimes, hence = 0. On mathematical ground one can argue that in the Hamiltonian 
framework, for the stationary black hole solutions, the total Hamiltonian function (which generates evolution 
along £ M ), must vanish as a first class constraint on the phase space [8,9]. This gives = E ADM . This is 
exactly what has been proved in another manner in the literature : that in stationary black hole spacetimes 
the ADM mass equals the energy of the black hole. Hence it is logical to identify with the horizon mass 
M b in the stationary case. The difference for an arbitrary non-stationary case is that H* ^ 0. Thus it can be 
called as the mass associated with the Isolated Horizon in an active sense, that can change from one dynamic 
equilibrium situation to another satisfying the first law. Here, one should be careful that this mass associated 
with the isolated horizon is completely physical and is not to be confused with the Hamiltonian of the SU(2) 
Chcrn-Simons theory which describes horizon field degrees of freedom. The Hamiltonian for Chcrn Simons 
theory vanishes identically, since the theory is topological and insensitive to small metric deformations. 

Clearly, the horizon mass is not affected by boundary conditions at asymptopia. It is defined locally on 
the horizon without referring to the asymptotic structure at all. The asymptotic conditions only modify the 
energy associated with the boundary at infinity and the bulk equation of motion (Einstein equations) [9,10]. 
The Hamiltonian framework discussed above is equally applicable for asymptotically flat and AdS spacetimes. 

2.2 Quantum Geometry 

For a classical spacetime with boundary, boundary conditions determine the boundary degrees of freedom and 
their dynamics. For a quantum spacetime, on the other hand, fluctuations of the boundary degrees of freedom 
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have a 'life' of their own (see for instance ref. [1]). Consequently, the Hilbert space of a quantum spacetime with 
boundary has the tensor product structure H = 1-L V <8>%b, with the subscript v (b) denoting the bulk (boundary) 
component. Thus, any generic state in quantum geometry, |\&), admits the expansion 

|*) = C vb |^> 8> \Xb) ■ (2.2.1) 

v.b 

In presence of electromagnetic fields, one can consider \ijj v ) (resp. \Xb}) to be the composite quantum 
gravity + quantum electrodynamics bulk (resp. boundary) state. The bulk states are annihilated by the full 
bulk Hamiltonian : H v \ip v ) = [H g , v + H e , v ]\ip v ) = 0; this is the quantum version of the classical Hamiltonian 
constraint [28]. The total Hamiltonian operator acting on the generic state |^) has the form 

# T |*> = {Hv ®h + h® Hb)\^) (2.2.2) 

where, I v (Ib) corresponds to the identity operator on T~i v {}i b ). 

While defining the grand canonical partition function, the charge operator (Q) for the black hole is also 
needed. It can be written in a similar fashion like the Hamiltonian as 

Ql*) = (Qv ®i b + iv® Qb)\v) 

where Q v and Q b are corresponding charge operators for the bulk states \ip v ) and the boundary states \xb) i 
respectively. In the classical theory the charge of a black hole is defined on the horizon i.e the internal boundary 
of the four dimensional spacetime (e.g. one can see how charge can be properly defined for spacetimes admitting 
internal boundaries in Einstein-Maxwell or Einstein- Yang-Mills theories in [7]). There is no charge associated 
with the bulk black hole spacetime, i.e. Q v , which is basically the Gauss law constraint for electrodynamics. 
Hence, its quantum version is assumed to be 

Qv\1>v) = 0. 

Combining the quantum constraints on the Hamiltonian and charge operators we can define a new quantum 
constraint as 

where H' v = Htv — &Qv and $ may be any function. But in our case it is a physically significant quantity 
which will be defined in the next paragraph. The implications of these quantum constraints will be seen during 
the construction of the grand canonical partition function. 

2.3 The Partition Function 

Let us consider a grand canonical ensemble of massive charged black holes immersed in a heat bath at some 
finite temperature with which it can exchange energy and charge as well. We construct a partition function for 
the thermodynamic system. Using the usual definition of the grand canonical partition function we write 

Z G = Tr exp -/3H T + (3$Q 

where the trace is taken over all states. $ is the electrostatic potential and Q is the charge operator for the 
black hole. To write it in the explicit form first we write down a general quantum state of the black hole as 
follows 

Now, we can write the partition function as 

Z G = Y,\^b\ 2 (Xb\®(^v\exp-/3H T + f3$Q\'ilj v )®\xb) 

v,b 

= Yl M 2 (»>l®<^l exp -pH'\ilj v )®\xb) (2.3.1) 

r.b 

where, H' = Ht — 3>Q- Writing the new operator H' as (H' v ® J& + I v ® H' b ) and using H' v \ip v ) = 0, the 
partition function comes out to be equal to the boundary partition function only i.e. Zq = Zq^, where Zcb is 
the boundary partition function for the charged isolated horizon, given by 

Z Gb = Tr b exp-f3(H b - $Q b ) 
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where it is assumed that the boundary states can be normalized through the squared norm \ c vb\ 2 (^vlipv) = 
\Cb\ 2 ■ This is analogous to the canonical ensemble scenario described in [19]. 

Now, the spectrum of the boundary Hamiltonian operator is still unknown in Loop Quantum Gravity. So 
we assume that the spectrum of the boundary Hamiltonian operator is a function of the discrete area spectrum 
and the charge spectrum of the area and the charge operators associated with the horizon, respectively 1 . The 
charge spectrum is equispaced on general physics grounds of charge quantization. The area spectrum, in Loop 
Quantum Gravity, can be approximately taken to be equispaced for large area black holes. For large area 
black holes, it is a well known fact in literature [31] that the major contribution to the entropy comes from 
the lowermost spins. Of course the higher spins contribute, but their contribution is exponentially suppressed. 
Thus, the effect of the corrections due to the inclusion of these higher spins is much more a technical aspect 
rather than physical. Hence, only spin 1/2 contributions for all punctures is taken into account, which leads to 
the equispaced area spectrum as an approximation, linear in the number of punctures . 

In a basis in which both area and charge operators are diagonal, the partition function can be written as 

Z G = ^2g(k,l)exp-p[E(A k ,Qi)-^Qi] (2.3.2) 

k,l 

where g(k, I) is the degeneracy corresponding to the area eigenvalue Ak and the charge eigenvalue Qi. k,l are 
the area and charge quantum numbers respectively. As we are interested in regime of the large area and charge 
eigenvalues (J; > 1,/ > 1), Poisson resummation formula is applied [17] to approximate the summation to an 
integration given by 

Zg = Jdx dy exp -p {E [A(x), Q(y)} ~ $Q(y)} g [A(x), Q{y)] (2.3.3) 

where x and y are area and charge quantum numbers in the continuum limit of k and I respectively Since 
A = A(x) and Q = Q{y), we can write dx = and dy — ^ to write the partition function in terms of area 
and charge as free variables as follows 

Zg = J 9(A, Q) exp -p {E(A, Q) - <S>Q} 

« JdAdQ e 8W-mA,Q)+f>*Q (2.3.4) 

where = e S(A) [ 25 ]. S{A) being the microcanonical entropy is a function of black hole area alone, as has 

been established within loop quantum gravity [1,3,4]. 



3 Stability Against Gaussian Thermal Fluctuations 
3.1 Saddle Point Approximation^. P.A.) 

Having a well defined partition function, we investigate its fmitcness under Gaussian thermal fluctuations 
about stable equilibrium configurations of the black hole given by the saddle points {A, Q}. Taylor expanding 
(S(A) — (3E(A, Q)+P&Q) about a saddle point (A, Q) and applying the saddle point conditions one can rewrite 
the partition function as 

Zg = e [S(A)-f3M(A,Q)+/3<!>Q] 

-±[-{S AA {A)-pM AA (A,Q)}a 2 +l3M QQ (A,Q)q 2 +2l3M AQ (A,Q)aq} da ^ (3.1.1) 

where, we have used M(A, Q) to indicate the equilibrium isolated horizon mass as a function of the area and 
charge. The saddle point conditions imply that the coefficients of a — (A — A) and q = (Q — Q) must vanish, 
which yield 

P(A Q) = J^f Q) , HA Q) = M Q (A, Q) (3.1.2) 



1 Actually this second assumption follows from the discussion in Subsection(2.1) [7,9] for spacctimcs admitting weakly isolated 
horizons where there exists a mass function determined by the area and charge associated with the horizon. This is an extension 
of that assumption to the quantum domain. 
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3.2 Quantum Surface Gravity 

An interesting result of this statistical mechanical approach is that, it gives rise to a quantum correction to 
the surface gravity which is a direct consequence of the logarithmic corrections of the microcanonical entropy 
5 = 4 ~ | l°g 4 from loop quantum gravity. If one calculates (3 from the saddle point conditions and use it 
to find the quantum surface gravity (K qU antum) m terms of classical surface gravity (^classical), there appear 
additional correction terms. One does this by calculating j3 from (3.1.2) and then using it in the expression 
^quantum = (more appropriately, this /3 can be replaced by (3 qU antum )■ The quantum corrections to the 
classical surface gravity is found to be 

where higher order terms are neglected for large black holes {A >> 1), A being the area of the weakly isolated 
horizon in Planck units. One can easily check the formula by applying it to the Reissner-Nordstrom and AdS 
Reissner-Nordstrom cases. Since the formulation is not dependent on any particular situation (symmetry, etc.), 
it will also be valid for other massive charged black holes also. 

3.3 Stability Criteria 

For the integral (3.1.1) to be convergent, the Hessian matrix H, given by (3.3.1) has to be positive definite 
where 

_ f/3M AA (A, Q) - Saa(A) (3M AQ (A, Q)\ 

{ pM AQ {A,Q) (3Mqq(A, Q)J {6 - 6 - L) 

The necessary and sufficient condition for the real symmetric square matrix H to be positive definite can 
be stated as follows [29] - 'Determinants of all the principal submatrices, including the determinant, of H are 
positive. ' It is also a crucial point that the inverse temperature f3 must be positive. Hence, the necessary and 
sufficient conditions for the positive definitencss of the Hessian matrix lead to the following stability criteria 

PM AA {A, Q) - S AA {A) > (3.3.3) 
detH = {pM AA (A,Q)-S AA {A)}pM QQ {A,Q)-p 2 M 2 AQ {A,Q)>0 (3.3.4) 

It should be noted that, for Mqq(A, Q) > 0, it will suffice to check only conditions (3.3.2) and (3.3.4) (which 
will be the case for RN and AdS RN black holes). 

Here, one may wonder that how these stability criteria (3.3.2)-(3.3.4) are related to the convexity property 
of the entropy function, which is the usual notion for thermodynamic stability. It is true that the usual notion 
of thermodynamic stability is related to the convexity property of the entropy function. It is also true that 
this convexity property of the entropy function follows from the requirement of the convergence of the partition 
function under Gaussian thermal fluctuations [17,22,25]. Our stability criterion, likewise, follows from equations 
(3.3.2)-(3.3.4) which constitute the necessary and sufficient conditions for the grand canonical partition function 
to be well defined. Physically, the conditions to be satisfied by the entropy Hessian in [17], [22] and [25] (which 
imply convexity of entropy function) and the conditions (3.3.2)-(3.3.4) lead to the same conclusion i.e. the 
finitcness of the partition function under Gaussian thermal fluctuations. In fact, one can check that the apparent 
dissimilarity in the mathematical structure of the Hessian in [17] and eq. (3.3.1) is just a manifestation of the 
different variables which are summed over 2 . The equivalence between our conditions and the convexity of the 
entropy function is obvious for the charge-less case i.e. canonical ensemble, discussed in details in [19]. 

3.4 The Product Ansatz 

The stability criterion (3.3.4) derived from the positive definitencss of the Hessian matrix, takes the form 

M AA M$ A S AA 

M A M QQ M A S a { ■ ■ ' 



To understand the motivation behind this choice of variables one can see Subsection (2.1) where a detailed explanation is given. 
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This is a second order partial differential inequality; one can resort to solving it for M(A, Q) in terms of S(A) 
and Q by trying to separate variables as one does for second order partial differential equations. Accordingly, 
we decompose M(A,Q) = n(A) ■ x(Q) with dim /it = dimM and dim% = 0. Further, x(0) = 1, so as to agree 
with the result in the chargeless (Schwarzschild) case. The inequality (3.4.1) then reduces to 



JL 



H>AA _ SaA 

~~s7 



/'-l 



> 



XQ XQ 
X XQQ 



(3.4.2) 



Since the left and right sides of inequality (3.4.2) are functions of two completely independent variables A 

and Q, for the inequality to make sense, one may set the right side (say) to a constant l/n say. One then 
obtains 

X(Q) = {1 + CQ)^ (3.4.3) 
where, C is a constant with dimC = dimQ . Substituting this in the inequality (3.4.2), we get 

n(A) > (aS + 7) 7r ^ T (3.4.4) 
where, k > 1. Setting 7 = for simplicity, and choosing (^a) = Mp, we get 



M 

Up 



> 



kB 
S_ 
k B 



k B (l + CQ) 
1 



1 



K — 1 



In 



S 



k B (l + CQ) 



(3.4.5) 



where the second line follows from the first in the limit k 1. It is easy to see that this bound is consistent 
with the chargeless case. It has to be seen how this bound checks out when the classical mass spectrum as a 
function of the area and charge of RN and AdS-RN black holes are used. The mass bound (3.4.5) is rewritten 
as 



M(A,Q) > S(A) 



S(A) 



l/re-l 



(3.4.6) 



where we have replaced the constant C in terms of lc- The classical RN is given by the metric (G = ks = c = 
h = 1) 



ds 1 



(1 + ^dt 2 + (1 + -V) dr 2 + r 2 dQ 2 



(3.4.7) 



This leads to the mass associated with the RN horizon as a function of the horizon area and charge: 



M RN (A,Q) 



-in 



1/2 r 



4ttQ 2 
A 



(3.4.8) 



Clearly, the charge is bounded by above according to 47rQ 2 < A. Inserting this in equation (3.4.6) leads to 



M(A) > S 



1 + SV2 



I/k-1 



(3.4.9) 



where, we have used 5(A) = A/ Alp and set lc = 2lp for simplicity. We have reverted to the use of the Planck 
length lp = (Gh/c 3 ) 1 / 2 as a fundamental length scale signifying onset of LQG. For black holes of very large 



horizon area A >> lp, the above bound reduces to 



M{A) > 



(3.4.10) 



In the same limit the classical Mrn obeys the upper bound Mrn < 2S 1 ' 2 which obviously violates the 
lower bound (3.4.10). 
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A similar analysis for the AdS-RN shows that for black holes of very large area compared to the cosmic area 
Aa = 4-7T (— t) = 47rZ 2 , the classical black hole mass obeys the bound (3.4.10). This observation comes from 
the fact that the metric for AdS-RN is 

, o . 2M Q 2 A o. , 9 , 2M Q 2 A 9n , , , , - 

rfs 2 = - 1 + -V - T r ) d * + 1 + ^7 " T-r 2 )- 1 ^ 2 + rW 3.4.11 

r H 3 r r z 3 

where, f = The mass associated with the AdS RN horizon is given by 



M ADS rn(A,Q) = [ — 



A^V 2 



A 4ttZ 2 



(3.4.12) 



It is obvious that for A » Aa, the last term in (3.4.12) dominates over the other two, and in this case, the 
bound is indeed obeyed. In fact the charge term makes a small contribution in any event, being bounded from 
above by 1. 

First Law : The mass formula written in the form M(A, Q) = p(A)x{Q) satisfies the first law. Differenti- 
ating M we have 

dM = HAXdA + [ixqAQ 
Now, from (3.1.2) we have $ = Mq = fJ-XQ- Hence 

dM = [lAXdA + $dQ 

Again, finding (3 from (3.1.2) and using n sg = one finds that the surface gravity comes out to be n sg = 
87r/i^x('sg' stands for surface gravity). Then, it is easy to see that this n sg proves that the product ansatz 
satisfies the first law 

dM = ^-dA + $dg 

Thus, from the thermodynamic point of view, this kind of solution for mass is acceptable. 
3.5 The Classical Metrics 

Reissner-Nordstrom black hole : In this paragraph we investigate the stability of massive charged Reissner- 
Nordstrom black holes against Gaussian thermal fluctuations. The classical Reissner-Nordstrom metric which 
is given by (3.4.7). The mass of the black hole in terms of area and charge as independent variables and it is 
given by 

M=i(^) 5 (l + P) (3.5.1) 

where p = w jf is a dimensionless parameter in the chosen units. The microcanonical entropy including the 
logarithmic correction term from loop quantum gravity is given by 

S=T-h°4 (3-5-2) 



Using (3.5.1) and (3.5.2) one finds that 



4 2 ° 4 



l 3 = 2^A-i x — -4f (3.5.3) 
(1-PJ 

detH = _-(izil£±i) 0.5.4) 

2A (1-p) K ' 

where p = 4 *9 . Conditions (3.3.2) and (3.5.3) together imply p < 1, whereas, (3.3.3) and (3.5.4) imply p > 1. 
From this contradiction one can conclude that the Reissner-Nordstrom black hole is locally unstable against 
Gaussian thermal fluctuations. 3 



3 This result is in agreement with what has been told in [24] i.e. an asymptotically flat, non- extremal black hole can never 
achieve a state of thermal equilibrium. This paper also uses a statistical mechanical formulation without any classical metric but 
do not have a sound mathematical basis as a justification of the quantum statistical formulation which has been presented in our 
work. 



7 



AdS Reissner-Nordstrom black hole : In this paragraph we investigate the stability of massive charged 
AdS Reissner-Nordstrom black holes against thermal fluctuations. The classical AdS Reissner-Nordstrom metric 
which is given by (3.4.11). The mass of the AdS Reissner-Nordstrom black hole in terms of area and charge as 
the independent variables is given by 



M 



-) 

4tt J 



(1 + P + *) 



(3.5.5) 



where we have set — = — jz as A is negative for AdS spacetimes and introduced the parameter a = j£p = . 
Here also the entropy is given by (3.5.2). Calculating the required quantities for AdS Reissner-Nordstrom case 
from (3.5.2) and (3.5.5) one finds that 



= 2^ A? 



(1-P+ 3a) 



detH 



2A (l-p + 3a) 2 



p-1 



1+T 



where a 



A 

4tt/ 2 



For the conditions (3.3.2) and (3.3.4) to be satisfied simultaneously the parameters have 



the following bound given by 



< p < 1 + 3(7 



which in the first order approximation reduces to the bound 



1 - 3ct 1 



24 

T 



< p < 1 + 3a 



Hence, AdS Reissner-Nordstrom black holes are locally stable against thermal fluctuations for the above range 
of parameters. 



4 Validity of S.P.A. 

To ensure the validity of the saddle point approximation, the relative r.m.s. fluctuations about the saddle point 
are to be checked i.e. AA rms /A = (AA 2 ) 1 / 2 / A, AQ rms /Q = (AQ 2 ) 1 ^ 2 /Q where 

i2_,u--i>, _ P M QQ \rfl-fTT-i\ _/3M aa -S A a 



AA = (H~ )u = ^-rrr , AQ 2 = (H-% 2 = 



detH v ' detH 

4.1 The Product Ansatz 

Using the product ansatz the mean square fluctuation of the charge comes out to be 

A q 2 = (SaPaa ~ S A aHa)paX 2 / 4 x ^ 

S a {SaPaa- S A aPa)hxqqX~ S 2 A p A x 2 Q 

Now, the following quantities are calculated 

C 2-„ $ 

XQ = + = _ 

K — 1 fJ, 



C 2 {2~k). -,^v»=2a C 2 (2-k) /k-1\^^3^ 
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Using all these quantities in (4.1.1) one finds that 



v/S/ 
77/1(2 - k) - S A n 2 j 



AQ 2 = —o 'T c.2 (4-1-3) 



where 77 = Sa/J-aa — Saa^a- Now, condition (3.3.4) implies Mqq((3Maa — Saa) > (3M A q. Since the R.H.S. is 
positive definite (/3 > 0), the L.H.S. has to be positive definite. This implies either Mqq > , ((3Maa — Saa) > 
or Mqq < , (/3Maa — Saa) < 0. The second one is not possible as it makes AQ 2 negative. Hence the only 
possibility is the first one. Since Sa > 0, hence ((3Maa — Saa) > implies 77 > 0. Keeping all these conditions 
in mind if one goes back to the expression (4.1.3), one finds that for AQ 2 > 0, k must satisfy the condition 
k < 2. The relative r.m.s. fluctuation of the charge goes as 



(AQ) rms 4C7T3 + M 



{(1 + | )(2 _ K) + (1 _6 )}3 ^ 



A~* (4.1.4) 



where /1 and its derivatives are calculated from the relation A = 8tt^, 2 . The arguement is as follows. The 
irreducible mass is a function of area only and it is the Schwarzschild mass. Hence, for Q = 0(x(0) = 1) we 
have A = 8nfi 2 . The mean square of area fluctuations comes out to be 

AA 2 = - ^ A 



TlfJl - 

Thus the relative rms fluctuation of area comes out to be 

(AA) 

rms I n^A 



Sa^ a 

2—K 



A In/; - Sa ^ / 167T/J 2 

To see the variation of (4.1.5) with A it is written explicitly in terms of A as follows 

{AA) rms 2V2 



(4.1.5) 



{(i-i) ^-Mi+Dr 



A--- (4.1.6) 



This again shows that the constant k has an upper bound which is not much larger than unity. Moreover, 
calculation yields 77 = — 64 ^- (l + "f)^ _5 which is a negative definite quantity(contradicts 77 > 0). For 77 < 
we must have XQQ < which in turn implies n > 2. But this step i.e. 77 < will lead to negative AQ 2 , which 
again implies instability. Hence k has to be bounded from above by 2. Considering the appropriate bounds 
of n it is seen from (4.1.4) and (4.1.6) that the relative r.m.s. fluctuations fall off for large area of the weakly 
isolated horizon i.e. A ^> 1 in Planck units. 



4.2 The Classical Metrics 

Reissner-Nordstrom black hole : The mean square fluctuations of area and charge in this case come out 
to be negative due to the presence of the determinant of the Hessian matrix. This emphasizes the instability of 
Reissner-Nordstrom black holes against Gaussian thermal fluctuations. 

AdS Reissner-Nordstrom black hole : The relative r.m.s. fluctuations of area and charge in this case 
arc found to be 



AA r 



A 



AQ r 



8(l-p + 3a) 



;i 

$2 



"*) {P- 
(l-p + 3a)(l-f){ 



A~ 



-(1-f) 



(4.2.1) 
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ir 1 } 



;i-d 



3-(l"f) 



I)" 1 } 



A~ 



(4.2.2) 
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where the equilibrium charge has been replaced in terms of the electric potential $ which can be easily calculated 
from the saddle point conditions (3.1.2). It is clearly seen from the above expressions that the fluctuations fall 
off in the large area regime. As far as the positivity of the mean square fluctuations(r.m.s. fluctuation to be 
real) is concerned one needs to be more careful with the lower bound of p. One can show that there arise the 
following two cases : 

1) For a < i (l + 21), the correct lower bound for p is what has been obtained i.e. 1 — 3<r (l — 21). 

2) For <7 > 1 (l + 21), the correct lower bound for p is given by 1 — 3<r (l — ^) + 

5 Discussion 

The first thing to say about this approach is that its origin is purely based on quantum aspects of spacetime. 
During the build up of the formalism, no classical metric is used. The construction of the partition function 
is purely based on the ideas and results of Loop Quantum Gravity e.g. the use of Chern-Simons states, the 
splitting up of the total Hilbert space, etc. and also on the Hamiltonian formulation of spacetimes admitting 
weakly isolated horizons. The entropy correction also follows from the quantum theory. The classical metrics 
come to the picture only to be tested. 

In course of this heuristic statistical mechanical approach of stability analysis of black holes, broadly two 
assumptions are made. In classical Hamiltonian GR it is known that the total Hamiltonian (gravity + matter) 
vanishes. So, it is very logical to consider that the quantum total Hamiltonian operator annihilates the bulk 
states of quantum matter coupled spacetime. A similar argument follows for the assumption of the quantum 
constraint on the volume charge operator. These two assumptions may be considered to be one due to their 
fundamental similarity and they ultimately give rise to a single quantum constraint. 

In Section(2), a second assumption is made regarding the eigenvalue spectrum of the energy of the black 
hole. It is already mentioned (in a footnote) that this second assumption is not a strong one. If one studies [7] 
carefully, the classical mass function associated with the horizon is stipulated to be a function of horizon area 
and charge. Again, this horizon area and charge are the functions of the local fields on the horizon. Proper 
quantization of the classical horizon area and charge will obviously lead to a well defined boundary Hamiltonian 
operator. The fact that there exists a quantum boundary Hamiltonian operator which acts on the boundary 
Hilbert space of the black hole is an assumption, since the exact form of such a Hamiltonian operator is still 
unknown. But the fact that its eigenvalue spectrum is a function of eigenvalue spectra of the area and charge 
operators is most likely a valid assumption, as it is bound to happen if such a boundary Hamiltonian operator 
exists. It follows from the classical analog - the mass associated with the horizon must be a function of the 
horizon area and charge for a consistent Hamiltonian evolution [7]. 

In ref. [17] where a similar approach has been taken, a particular functional form of the mass in terms of 
the area and charge had been used on an ad hoc basis. Such an ad hoc assumption has been shown here to 
be quite redundant. This, therefore, is a significant strength of the paper, relative to the earlier assay. Thus, 
the statistical mechanical approach adopted in this paper, though similar, now stands on a far stronger ground 
than in the previous version. 

This statistical mechanical approach gives us a new quantum correction to the surface gravity arising from the 
loop quantum gravity corrections to the microcanonical entropy. One can easily check its validity. Moreover, 
it predicts local thermodynamic instability of the Reissner-Nordstrom black hole and local thermodynamic 
stability of the AdS Reissner-Nordstrom black hole as was shown classically in [21]. As far as the relative charge 
fluctuations are concerned, in literature [23] there are problems for AdS RN black holes as it does not fall off 
for any condition. This problem has been solved in this papcr(Scction4). Last but not the least, once more a 
word is worth mentioning that this formulation does not involve any classical metric. The whole thing depends 
on quantum aspects of spacetime. One can generalize this to study rotating and even charged-rotating black 
holes because there is no use of symmetry in the theory. We look forward to give those analyses in future. 

There is a more crucial issue which can be of utmost importance to have a deeper understanding of black 
hole thermodynamics [30]. As far as the saddle-point approximation is concerned, the Euclidean path integral 
approach does look similar to our thermal holographic approach. But there is a crucial difference between the 
two approaches. The path integral approach, within the saddle point approximation, requires information of 
the full black hole spacetime (bulk and horizon) as given by the classical black hole metric solution chosen to be 
the saddle point. Thus one needs global information away from horizon. In the thermal holographic approach 
which we adopt, one needs only local information associated with the equilibrium isolated horizon geometry 
interpreted as an inner null boundary. Thus, no detailed knowledge of the full classical spacetime is needed. 
The mass of the isolated horizon is an unspecified function of the area and charge, and we never need to specify 



fl) 



this function to derive our results, except as fiducial checks [Subsection (3.5)] appropriate to given classical 
metrics. This insensitivity of our approach to an explicit classical black hole metric is a key feature of our work 
and can be taken to mean that our results are in a sense more general than those computed from the Euclidean 
path integral. 

A further distinction is that, in contrast to the Euclidean path integral approach, where quantum fluctuations 
around a classical metric are considered, our saddle point is a quantum isolated horizon whose quantum states 
and their (non-perturbative) dynamics are described by a quantum Chern Simons theory. Consequently, the 
equilibrium entropy is the microcanonical entropy computed in earlier work (ref. [1-4]) based on Loop Quantum 
Gravity, and already has an infinite series of corrections (including those logarithmic in horizon area) beyond the 
Bckcnstcin-Hawing area law, incorporating quantum spacetime fluctuations. In this paper, additional thermal 
fluctuations (and their physical effects) arc considered, over and above the quantum spacetime fluctuations 
already incorporated for the equilibrium configuration. Quantum and thermal fluctuations arc thus, treated 
somewhat distinctly in our approach, and the result is an interesting interplay between them. In the Euclidean 
path integral, such distinctions are not as clear. It might be of future interest to see better how these two 
somewhat disparate approaches may be related. 

Acknowledgments : We thank T. Sarkar for many illuminating discussions during the course of this work 
and also for bringing reference [22] to our attention. 
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